Chaos edges of ^-logistic maps: Connection between the relaxation and sensitivity 

entropic indices 
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Chaos thresholds of the z-logistic maps xt+i = 1 — a|a;t| z (z > l;t = 0,1,2,...) are numer- 
ically analysed at accumulation points of cycles 2, 3 and 5. We verify that the nonextensive q- 
generalization of a Pesin-like identity is preserved through averaging over the entire phase space. 
More precisely, we computationally verify ]im.t->ao(S q ™ n )(t)/t — Hmt->oo (ln g jj n £}(t)/i = X q av , 
where the entropy S q = (1 — ^2 i pj) / (q — 1) (Si = — ^2 i Pilapi), the sensitivity to the initial 
conditions £ = limA:ri(o)-> o Aa;(t)/Aa;(0), and ln 9 x = (a; 1-9 — 11/(1 — q) (hii x = \nx). The en- 
tropic index q^ n < 1, and the coefficient X q av > depend on both z and the cycle. We also 
study the relaxation that occurs if we start with an ensemble of initial conditions homogeneously 
occupying the entire phase space. The associated Lebesgue measure asymptotically decreases as 
l/tVCarei- 1 ) (q rel > 1). These results led to (i) the first illustration of the connection (conjectured 
by one of us) between sensitivity and relaxation entropic indices, namely q re i — 1 ~ ^4.(1 — qge„) a , 
where the positive numbers (^4, a) depend on the cycle; (ii) an unexpected and new scaling, namely 
qsLicycle n) = 2.5 q™ n {cycle 2) + e (e = -0.03 for n = 3, and e = 0.03 for n = 5). 
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Boltzmann-Gibbs (BG) entropy and corresponding 
statistical mechanics generically require strong chaos for 
their applicability and (notorious) usefulness. This type 
of requirement was first used in 1872 by Boltzmann him- 
self [1]. Indeed, his "molecular chaos hypothesis" al- 
lowed him to arrive to the celebrated distribution of en- 
ergies at thermal equilibrium, now known as Boltzmann 
weight. Today, we know that this requirement essentially 
amounts, for classical nonlinear dynamical systems, to 
having at least one positive Lyapunov exponent. In the 
case of many-body Hamiltonian systems, such a condi- 
tion is satisfied when the interactions are short ranged. 
Such systems typically exhibit three basic exponential 
functions [2], namely (i) the sensitivity to the initial 
conditions diverges exponentially with time, (ii) physical 
quantities exponentially relax with time to their value 
at the stationary state (thermal equilibrium), and (iii) 
at thermal equilibrium, the probability of a given mi- 
crostate exponentially decays with the energy of the mi- 
crostate. These three exponentials of different, though 
related, nature can be summarized in the following dif- 
ferential equation: 



dx 



ai y (y(0) = l) 



(1) 



whose solution is y = e aiX (the subindex 1 will become 
transparent soon). Let us make explicit the point. The 
first physical interpretation concerns the sensitivity to 
the initial conditions of say a one-dimensional case and 
is defined as 



Ax(t) 
Ax(ay->o Ax(0) ' 



where Ax(t) is the distance, in phase space, between 
two copies at time t. If the system has a positive Lya- 
punov exponent Ai, then £ diverges as £ = e Alt . In 
other words, in this case we have (x, y, ai) = (i, £, Ai). 
The second physical interpretation concerns the relax- 
ation of some (characteristic) physical quantity 0(t) to 
its value O(oo) at thermal equilibrium. With the defini- 
tion Q = p[o)Zp(^) , w e typically have Q = e _t / T , where 
r is the relaxation time. In other words, in this case 
we have (x, y, a±) (t, f2, — 1/r). This relaxation occurs 
precisely because of the sensitivity to initial conditions, 
which guarantees strong chaos. It was apparently Krylov 
the first to realize [3], over half a century ago, this deep 
connection. The third physical interpretation is given by 



Pi 



2^=1 1 



-0E 3 



), where Ei is the 



eigenvalue of the i-th quantum state of the Hamiltonian 
(with its associated boundary conditions), pi is the prob- 
ability of occurrence of the i-th state when the system is 
in equilibrium with a thermostat whose temperature is 
T = l/k/3 (Gibbs' canonical ensemble). In other words, 
in this case we have (x, y, a\) — > (Ei, Zpi, —j3). 

A substantially different situation occurs when the 
maximal Lyapunov exponent vanishes. In this case the 
typical differential equation becomes 



dy 
dx 



.y* (y(Q) = l;qeK), 



(3) 



whose solution is y = e a q ' 
being defined as follows: 



the q-exponential function 
= [1 + (1 - g)x] 1 /( 1 -?) if 



lim 



(2) 



the quantity between brackets is nonnegative, and zero 
otherwise (ef = e x ). The sensitivity to the initial condi- 



tions is given in this case by [4,5] £ = 



(sen stands 



1 



for sensitivity). In other words, we have (x,y,q,a q ) = 
{t,£,q se n,\ aen )- The relaxation is typically expected 

—t/r 

[6] to be characterized by = e qrlll " r ° l (rel stands 
for relaxation). In other words, in this case we have 
(x,y,q,a q ) = {t,Cl,q re i,-l/r qrel ). For the longstanding 
metastable states [7] that precede thermal equilibrium 
for long-range interacting Hamiltonians, it is expected 

[8] Pl = e-f t :r tEi /z gstat (with z qstat ^ yJ=i e q f t :r tEi ) 

{stat stands for stationary). In other words, in this case 
we have (x,y,q,a q ) = (Ei,Z qstatP i,q stat ,-f3 qstat ). For 
systems that are at, or close to, the edge of chaos we typ- 
ically have q sen < 1, q re i > 1, and q stat > 1. For the BG 
case, where there is one or more positive Lyapunov expo- 
nents, we recover the confluence q sen = q re i = q s tat — 1- 

One expects the entire g-triplet to be either measurable 
or calculable for Hamiltonian systems. And indeed it has 
recently been measured in the solar wind [10]. However, 
the generic relation among these three q indices is still un- 
known. For dissipative systems such as say the z-logistic 
map, no q s t a t exists. Therefore, the problem reduces to 
only two q indices, namely q sen and q re i. Their generic 
relation also is unknown. In the present paper, we pro- 
vide the first (numerical) evidence of such a connection. 

Before entering into the details of the present calcula- 
tion, let us briefly review the connection with the entropy 
S q , basis of a current generalization of BG statistical me- 
chanics referred to as nonextensive statistical mechanics 
[11]. This entropy is defined as follows: 



S q = 



1 



^i^ = 5>ln 9 (lM) 



(4) 



exponential, is defined as ln q x = x 1 q _ q 1 (hii x = lnx) 

and Si = Sbg = - YnLiPi ln Pi- 

If we partition the phase space of a one-dimensional 
map (at its edge of chaos) into W small intervals, ran- 
domly place N initial conditions into one of those win- 
dows, and then run the dynamics for each of those N 
points, we get, as time t evolves, an occupancy char- 
acterized by {Ni(t)} {EZi N i(t) = N )- With Pi(*) = 
Ni (t) /N we can calculate S q (t) for any value of q. From 
this, we can calculate the entropy production per unit 
time [12], defined as follows: 



lim lim lim 

t— >oo W— >oo AT— >oo 



s g (t) 



(5) 



where the q-logarithm function, inverse of the q- 



It has been proved [13] that only K qeen is finite (K q = 
for q > q sen and K q diverges for q < q sen ). Further- 
more, if we consider the upper bound of K qgcn with re- 
gard to the choice of the little window within which we 
put the N initial conditions, we obtain the Pesin-likc 
identity K qst . n — \ qscn . Several aspects of this problem 
have already been verified for various one-dimensional 
unimodal maps [14-17]. It was recently studied the in- 
fluence of averaging [18]. It was verified that, while the 
g-generalized Pesin-like identity is preserved, the value 
of q sen is changed into q™ n (av stands for average). The 
main goal of the present paper is to exhibit that a simple 
relation exists between q^ n and q re i by making use of 
the 2-logistic map family defined as 

x t +i = 1 - a\x t \ z (6) 

where (z > 1; < a < 2; \x t | < 1; t = 0, 1, 2, ...). 



TABLE I. z-logistic map family for cycles 2, 3 and 5. 
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cycle 




av 
Hsen 


qrei 




1&£2L 


1.75 


2 


1.355060... 


0.37 ±0.01 


2.25 ± 0.02 


0.26 ±0.01 


0.26 ±0.02 


1.75 


3 


1.747303... 


0.92 ±0.01 


2.25 ± 0.02 


0.48 ±0.01 


0.47 ±0.02 


1.75 


5 


1.607497... 


0.96 ±0.01 


2.25 ± 0.02 


0.42 ±0.01 


0.40 ±0.02 


2 


2 


1.401155... 


0.36 ±0.01 


2.41 ± 0.02 


0.27 ±0.01 


0.27 ±0.02 


2 


3 


1.779818... 


0.88 ±0.01 


2.41 ± 0.02 


0.49 ±0.01 


0.48 ±0.02 


2 


5 


1.631019... 


0.93 ±0.01 


2.41 ± 0.02 


0.42 ±0.01 


0.40 ±0.02 


2.5 


2 


1.470550... 


0.34 ±0.01 


2.70 ± 0.02 


0.28 ±0.01 


0.28 ±0.02 


2.5 


3 


1.828863... 


0.82 ±0.01 


2.70 ± 0.02 


0.48 ±0.01 


0.47 ±0.01 


2.5 


5 


1.669543... 


0.88 ±0.01 


2.70 ± 0.02 


0.38 ±0.01 


0.37 ±0.01 


3 


2 


1.521878... 


0.32 ±0.01 


2.94 ±0.02 


0.29 ± 0.02 


0.29 ±0.03 


3 


3 


1.862996... 


0.78 ±0.01 


2.94 ±0.02 


0.44 ±0.01 


0.44 ±0.01 


3 


5 


1.699440... 


0.84 ±0.01 


2.94 ±0.02 


0.34 ±0.01 


0.35 ±0.01 


5 


2 


1.645533... 


0.28 ±0.01 


3.53 ± 0.03 


0.30 ± 0.02 


0.30 ±0.03 


5 


3 


1.931072... 


0.68 ±0.01 


3.53 ± 0.03 


0.36 ±0.01 


0.37 ±0.01 


5 


5 


1.773088... 


0.73 ±0.01 


3.53 ± 0.03 


0.27 ±0.01 


0.25 ±0.02 
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FIG. 1. The volume occupied by the ensemble as a func- 
tion of discrete time. After a transient period, which is the 
same for all N box values, the power-law behavior is evident. 
For each case, the evolution of a set of WNbox identical copies 
of the system is followed. 

In our simulations, we firstly check whether q re i val- 
ues of cyles 3 and 5 are different from those of cycle 2 
obtained in [6]. To accomplish this task, we analyse the 
rate of convergence to the critical attractor when an en- 
semble of initial conditions is uniformly distributed over 
the entire phase space (the phase space is partioned N box 
cells of equal size) and we found that, for all cycles that 
we studied, the volume W(t) occupied by the ensemble 
exhibits a power-law decay with the same exponent value 
for fixed z. As an example, the case of cycle 3 for z = 2 
is given in Fig. 1. The same kind of behavior is ob- 
tained also for other z values and cycles, which yields us 
to conclude that q re i does not depend on the cycle (see 
also the Table). Then, we concentrate on the ensem- 
ble averages of the sensitivity function £(£) by consider- 
ing two very close points (throughout this work we take 
Ax(0) — 10~ 12 ) and calculating its value from Eq. (2). 
This procedure has to be repeated many times with dif- 
ferent x values randomly chosen in the available phase 
space and finally an average is taken over all ln g £ (t) val- 
ues. For cycle 3 and cycle 5, we obtain the behavior of 
(ln 9 £)(i) as a function of t, for various z values, from 
where one can deduce q^ n by identifying the linear time 
dependence as it is seen in Fig. 2. We verify that q™ n and 
A^a„ values do depend on both z and the cycle, whereas 
q re i is independent of the cycles. Finally, to investigate 
the entropy production for cycles 3 and 5, we employ the 
procedure used so far in [16] for cycle 2 of the z-logistic 
map family. It is numerically verified that, as seen for 
a representative case in Fig. 3, for each z value and for 
each cycle, the linear entropy production occurs for a spe- 
cial q value which coincides with the one obtained from 



the sensitivity function. In addition to this, we obtained 
= A°a„ , which clearly broadens the validity region 
of the standard Pesin theorem [19]. 
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FIG. 2. The behavior of {ln q £) as a function of time. 
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FIG. 3. The behavior of {S q ) as a function of time. 

Finally, for all cycles that we studied, we numerically 
verified that a simple scaling relation exists between q 
and q re i (see Fig. 4): 



av 
sen 



q re l(cycle n) - 1 ~ A [1 



q a J n (cyden)] c 



(7) 



where n = 2, 3, 5 and the values of A and a are given in 
the caption of Fig. 4; both numbers depend on the cycle. 
For example, a = 5.1 for cycle 2, and quickly approaches 
zero when the cycle increases; A also decreases when the 



3 



cycle increases. This kind of relation between these two 
classes of q index is seen for the first time in a model sys- 
tem. It is clearly consistent with the confluence occurring 
for BG systems. This is to say, when there is at least one 
positive Lyapunov exponent, we obtain q re i — q^ n = 1. 

We also notice (see Fig. 5) a new and unexpected scal- 
ing behavior, namely 

CJcyde n) = 2.5 q™ n (cycle 2) + e , (8) 

with e = —0.03 for n — 3, and e = 0.03 for n = 5. 

Summarizing, we have discussed a paradigmatic fam- 
ily of one-dimensional dissipative maps, and have shown 
that its (averaged) sensitivity to the initial conditions 
and its relaxation in phase space follow a simple path, 
which is consistent with current nonextensive statistical 
mechanical concepts, and which considerably extends the 
validity of Pesin-like identities. The sensitivity to the ini- 
tial conditions is characterized by q^ n < 1, which mono- 
tonically approaches unity with increasing cycle size (at 
least for the specific cycles that we have studied here), 
and decreases with z. It is further characterized by \®% v , 
which exhibits a maximum both as a function of the cycle 
size and of z. The relaxation is characterized by q re i > 1, 
which monotonically increases with z and does not de- 
pend on the cycle. This study has enabled to exhibit 
two interesting relations, namely Eqs. (7) and (8). This 
path is expected to appreciably enlighten, among others, 
the case of long-range- interacting Hamiltonian systems, 
where the situation is even more complex since a third 
entropic index, q s tat, is expected, which would character- 
ize the energy distribution at metastable states. Analytic 
analysis of the scalings presented here are certainly most 
welcome. 
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<T sen(cyc ie2) 

FIG. 5. Straight lines: qsen(cycleS) = 2.5 q^n {cycle 2) 
- 0.03, and q a J n (cycle 5) = 2.5 q™ n {cycle 2) +0.03, which sug- 
gests q™ n (cycle5) - qsen{cycle3) ~ 0.06. 
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FIG. 4. Straight lines: q re i(cycle2) 

? n {cycle2)] 5 - 1 , q rel (cycle3) - 1 = 4.6 [1 - q^ n (cycle3)f 
q rel (cycle5) - 1 = 4.1 [1 - g s 7 n («,de5)] - 39 . 
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